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We have studied the Junction Magnetoresistance (JMR) and the Differential junction magnetoresis- 
tance (DJMR) for double tunnel junctions with magnetic metals in the Coulomb Blockade regime. 
Spikes are seen in both the JMR and the DJMR vs. voltage curves. They occur at those places 
where the current increases by a step. In all cases the large bias limit can be obtained by adding 
the resistances of each of the junctions in series. The JMR is positive in all the cases we studied, 
whereas the DJMR can be positive or negative as a function of the voltage. Moreover, the relative 
variation of the DJMR as a function of the voltage is larger than the variation of the JMR with the 
voltage. 

PACS numbers: 73.40. Gk, 73.23.Hk, 75. 70. Pa 



I. INTRODUCTION 

Electron tunneling through tunnel junctions which 
consist of magrietic metals separated by insulators has 
a long historyJd Some of the earliest experiments were o: 
tunneling between a superconductor and a ferromagnet. 
In these experiments, the difference between the spin-up 
and spin-down density of states of a superconductor in 
a magnetic field was observed. These experiments also 
provided a quantitative measure of the spin polarization 
of the current from the ferromagnetic metal. Later ex- 
periments were done on tunneling between two ferromag- 
netic metals.0 Because of the spin dependence of the den- 
sity of states, the conductance changed when the relative 
orientation of the two magnetizations changed. Tunnel- 
ing between two ferromagnetic metals has experienced 
renewed interest because of its potential application in 
information storage. One of the measurable quantities of 
practical importance is the junction magnetoresistance, 
which is defined as the ratio of the change in the con- 
ductance from parallel to antiparallel alignment divided 
by the parallel conductance. Large magnetoresistance 
response in low magnetic fields has been obaepted in a 
variety of tunnel junctions by several groups.oilil 

Besides spin, electrons carry charge and at low tem- 
peratures in tunnel junctions with small capacitance, the 
discreteness of electron charge manifes±s_itself through 
the physics of the Coulomb Blockade.EjEJ The transfer 
of an electron by tunneling between two initially neu- 
tral regions of capacitance C increases the electrostatic 
energy of the system by an amount e^/2C. At low tem- 
peratures and small voltages, the tunneling current is 
suppressed because of the charging energy. This energy 
barrier is called the Coulomb Blockade. The effect of 
this in a two-junction system is the incremental increase 
in the current at voltages where it is energetically favor- 
able for an electron to tunnel to the center island. The 
occurrence of these current steps in the current-voltage 



characteristics is known as the "Coulomb Staircase." Ad- 
vances in technology have led to a number of new experi- 
ments which exhibi^Ji^ Coulomb Blockade phenomenon: 



granular thin-fiimsElrnd, lithographically patternedJ^iu, 



nel junctionsEijO, narrow disordered quantum wire 
scanning tunneling microscopy of small metal dropletsEH, 
metal-insulator- metal tunnel junctions pSstith small metal 
particles embedded in the insulator .EJtll 

Recently there have been experiments where spin 
polarized tunneli*^ and the Coulomb Blockade have 
been combinedHa For example, in discontinuous 
metal/insulator multilayers consisting of closely spaced 
ferromagnetic nanoparticles in an insulating matrix, a 
Coulomb Blockade is seen in the current-voltage charac- 
teristic at low temperatures. The purpose of this paper 
is to further study the nonlinear conductances for paral- 
lel and antiparallel magnetization alignment for a double 
junction with magnetic metals in the Coulomb Blockade 
regime. 

The organization of the rest of the paper is as follows. 
In the next section we describe the model and the mas- 
ter equation used. The expression for the current, the 
conductance, the differential conductance and the junc- 
tion magnetoresistance are then introduced. The results 
of the calculation are given in the discussion of Sec. Ill, 
and summarized in Sec. IV. 



II. FORMALISM 

The experimental realization of a double tunnel junc- 
tion consists of three metals separated by oxide layers 
which act as insulators. The two end leads, denoted by 
left (i) and right (i?), are connected to a constant volt- 
age source V . This double junction can be modeled as 
a closed loop electrical circuit with two leaky capacitors 
of capacitances Cl and Cr, connected in series (Fig. 1). 
The resistances of the junctions, Rl and Ru, are con- 
nected in parallel to Cl and Cr respectively. Depending 
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FIG. 1. Schematic diagram of a double tunnel junction sys- 
tem. The left and right capacitors {Cl and Cr) are connected 
in series with a voltage source. The resistances of the left and 
right junctions are Rl and Rr, respectively. The voltage in 
the left and right leads are kept fixed at V and 0, whereas in 
the central region the voltage Vm changes. 

on the number of excess electrons, n, in the center metal, 
the voltage in the central region Vm fluctuates. The volt- 
age drops, (Vl — Vjvf (n)) and (VA/(f^) — Vr), across the 
left and the right junctions are 
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where — e is the charge of the electron, C is the capaci- 
tance equal to (Cl + CpXy-Mui Vg plays the role of the 
gate or the offset voltage. c^Eil To describe. lIi ejbransport 
of electrons we use a semiclassical model.EjOL3 The dif- 
ferent tunneling rates that enter into the dynamics arc: 
the rate for electrons to tunnel onto the center metal 
from the left (F^^^^^^) and right (r^i^„_|_j), and the rate 
for the electrons to leave the center metal to the left 
ij'n^n-i) ^^^d right (r^_^„_2)- These ijjnnel rates can 
be calculated using Fermi's Golden rule.Ej To express the 
different tunnel rates in a convenient fashion we introduce 
a function 7(e): 



7(e) = 



1 



(3) 



In terms of the above function the various tunnel rates 
are: 
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Here Eq = jlC is the charging energy, and /3 is related 
to the temperature T via /3 = l/ksT. 

For magnetic materials, there are two types of carriers: 
the majority electrons with spins parallel to the magneti- 
zation and the minority electrons with spins antiparallel 
to the magnetization. We make the assumption that the 
rate of tunneling through the junctions is much smaller 
than the rate of spin relaxation. Under this assumption 
the number of majority and minority electrons in the 
center metal is not conserved separately, but their sum 
is conserved. This results in a master equation which de- 
scribes the time evolution of the probability, Pn{t), for n 
excess electrons in the center metallic 



dt 



— [Pn+lFn+l^n ~ Pnr„_+„_|_i] 



(5) 



where the net tunnel rate Ti^j is equal to (Ffl,^ -f- Ff^^ ) . 
At steady state, dpn{t)/dt = 0, the current in the left 
junction equals the current in the right junction due to 
current conservation, and is given by 
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To evaluate the current in general, Eqs. (||) - (||) have 
to be solved numerically. We follow the same procedure 
described in Ref. [35]. At low temperature and voltage, 
the higher charge states |n| >> 1 are energetically for- 
bidden and therefore we can neglect them. We found 
that for our choice of parameters it is sufficient to keep 
the 19 states around n — 0. The highest and the lowest 
states are n = 9 and n — —9. The total conductance 
G is then simply given by dividing the current with the 
voltage, G = I/V, and the differential conductance Gdiff 
is obtained by differentiating the current with respect to 
the voltage, Gdiff = dl/dV. 

There are two possible choice of materials which lead 
to non-zero magnetoresistance in a two junction system: 
(1) all the leads are ferromagnets, i.e. FM/I/FM/I/FM, 
and (2) the left or the right lead is a paramagnet, i.e. 
FM/I/FM/I/PM or PM/I/FM/I/FM. Note that with 
the center lead non-magnetic (FM/I/PM/I/FM) there 
is no JMR or DJMR because the antiparallel conduc- 
tance is identical to the parallel conductance within our 
assumption. For each of these cases, we calculate the 
conductances for two different relative orientations of 
the ferromagnet magnetizations: parallel and antiparal- 
lel. Here parallel (P) means that the magnetizations of 
all ferromagnets is in the same direction, and antipar- 
allel (AP) means the magnetization of the center ferro- 
magnetic metal has been reversed. In case (1) of three 
ferromagnetic metals, other configurations are also pos- 
sible: the magnetization of the right or the left lead can 
be reversed keeping the magnetizations of the other two 
parallel. For case (2) where one of the left or the right 
lead is non-magnetic, the magnetization of right or the 
left ferromagnetic lead can also be reversed. However, 
the qualitative features of the current, conductance, and 
the magnetoresistance curves remain the same. The junc- 
tion magnetoresistance (JMR) is then defined as the ra- 
tio of the change in the conductance from parallel, G^ 
to antiparallel alignments, G^^, divided by the parallel 
conductance 
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The differential junction magnetoresistance (DJMR) is 
defined in the same fashion as the JMR but with the total 
conductances replaced by the differential conductances, 
Gdiff, 
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DJMR = 1 



(8) 



III. DISCUSSION 

The different parameters of our problem are: the ca- 
pacitances Cl and Cr^ the parallel resistances i?^ and 
i?^, the antiparallel resistances R^^ and R^ , the tem- 
perature r, and the gate or the offset voltage Vq- To esti- 
mate the antiparallel resistances, we use Julliere's model 
for FM/I/FM tunnel junctions, where the junction mag- 
netoresistance is expressed as a product of t]a,e magni- 
tudes of the polarizations of the ferromagnets 
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In the above equation, AP is the change in resistance 
from antiparallel to parallel orientations, R^^ is the junc- 
tion resistance when the magnetizations of the ferromag- 
nets are antiparallel, and Pi and P2 are the spin polar- 
izations of the two FMs. 

In Fig. (2) the currents and the total conductances 
are plotted for parallel and antiparallel orientations of 
the magnetizations. Asymmetric junctions and low tem- 
perature are used to obtain the steps in the current. Also, 
we take the offset or the gate voltage to be zero, Vq — Q- 
The effect of finite gate voltage is only to shift the I — V 
curves.E3 As expected in the Coulomb Blockade regime, 
the current increases by steps (case a), and the conduc- 
tance shows oscillations with decreasing amplitude (case 
b). Finally, for large enough voltage, the current becomes 
linear with voltage and the conductance approaches the 
classical limiting value which is the inverse of the resis- 
tance of two resistors connected in series 



G = 
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The effect of temperature dependence is shown in Fig. 
(^. With the rise in temperature, the number of acces- 
sible states increases. The steps in the current round off 
(case a), and the amplitude of the oscillations in the con- 
ductance vs. voltage curve (case b) decreases. The peaks 
in the JMR (case c) and the DJMR vs. voltage curves 
(case d) become broader and reduced in size. 

Restricting ourselves to ksT « Ec and zero gate 
or offset voltage {Vq = 0), in Fig. 4 we plot 
the junction magnetoresistance for four different cases: 
(a) FM1/I/FM1/I/FM2, (b) FM2/I/FM1/I/FM1, (c) 
FMi/I/FMi/I/PM, and (d) PM/I/FMi/I/FMi, where 
FMi and FM2 are different ferromagnetic metals. In the 
high bias regime for all fouf-cases, the current is linearly 
proportional to the voltages 
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FIG. 2. Current (I) and total conductance (G) as a function 
of voltage, for a FM/I/FM/I/FM tunnel junction. For both 
the cases (a) and (b), the solid lines represent the current and 
the conductance for parallel alignment of the magnetizations, 
whereas the dashed lines are for antiparallel alignment. As 
a function of voltage the current shows steps (case a), and 
the total conductance oscillates with decreasing amplitudes 
(case b). We choose the left and the center metals to be 
made of iron with polarization, Ppc=0.40, thejright metal to 
be made of cobalt with polarization Pco=0.34.H The different 
parameters used are: Cl=0.01C, -R£ = 
keT = O.OlPc- 



--0.01{Rl + Rl) 



and 



It then follows from Eq. (|ll|), that the plateaus in Figs. 
4(a) - 4(d) are given by the difference in the total resis- 
tances from antiparallel to parallel alignments divided by 
the total resistance for the antiparallel alignment. 
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For asymmetric junctions in the low voltage, low tem- 
perature regime we see spikes in the JMR vs. voltage 
curves in Fig. (4). The occurrence of the first spike can 
be understood analytically. We confine ourselves to only 
three states: the n = state and the two states n = ±1 
immediately accessible from it. At very low temperature, 
the state which is most likely occupied is the n = state 
for Vg = 0. The tunnel rates Fq^i and Fi^o do not enter 
in the expression for the tunneling current at zero tem- 
perature, but at finite temperature they do contribute. 
In the limit of zero temperature, where we replace the 
function 7(e) in Eq. (||) by 7(e) = e^(e), the steady state 
current can be expressed in terms of the tunnel rates in 
the regime (e/2) < CrV < (Cfl/CL)(e/2) as 
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^^^^ In the derivation of Eq. ( p^ ) we have assumed that the 
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FIG. 3. Temperature dependence of the (a) Cur- 
rent, (b) Conductances, (c) JMR, and (d) DJMR for a 
FM/I/FM/I/FM tunnel junction. As in Fig. 2(b), the solid 
and the dashed lines represent the current and the conduc- 
tance for parallel and antiparallel alignments. With increase 
in temperature, the steps in the current vs. voltage curve 
round off (case a), resulting in reduced amplitude of oscilla- 
tions in the total conductance (case b). The peaks in both the 
JMR (case c) and the DJMR (case d) curves broaden. Also 
the peaks in the JMR and the DJMR curves decreases with 
increase in temperature. The parameters are the same as used 
in Fig. 2, except for the temperature which is fcsT = O.l-Ec- 

capacitance of the right junction is larger than the ca- 
pacitance of the left junction, i.e, Cr > Cl- Near the 
special point CrV « e/2, where the current increases 
by a step at zero temperature, the height of the spike is 
simply given by the difference, 

JMR(C«1/=^,T-.0) = 1--^. (14) 

However, at finite temperature, the height of the spikes 
depend on the temperature, the capacitances, and the 
resistances of the junctions. We also notice that spikes 
go up for cases 4(a) and 4(c) whereas spikes go down 
for cases 4(b) and 4(d). Cases 4(a) - 4(b) or cases 4(c) 
- 4(d) differ only by the exchange of the metals in the 
left and the right leads. In all the above cases we keep 
the resistance of the left junction for the parallel align- 
ment fixed, i?^ ~ 0.01 (i?^ + Rr)- By only exchanging 
the metals of the left and the right leads, we go from 
spikes which go up to spikes which go down and vice 
versa. The interchange of the spikes from up to down or 
down to up can be explained in the following way. At 
zero temperature, whenever the JMR obtained from Eq. 
( p^ ) is greater than the height of the plateau (Eq. (p^), 
we see an up-spike (Fig. 4(a) and 4(c)), otherwise we 
see a down-spike (Figs. 4(b) and 4(d)). Within JuUiere's 



FIG. 4. Junction magnetoresistance (JMR) as a func- 
tion of voltage for different cases: (a) FM1/I/FM1/I/FM2, 
(b) FM2/I/FM1/I/FM1, (c) FMi/I/FMi/I/PM, and (d) 
PM/I/FMi/I/FMi, where FMi and FM2 are two different 
ferromagnets, and PM is any paramagnet. In all cases we see 
spikes in the JMR: spikes go up for cases (a) and (c) whereas 
spikes go down for cases (b) and (d). These spikes occur at 
those places where the current increases by a step (see Fig. 
2(a)). The plateaus in all cases are given by the difference in 
the total resistances from antiparallel to parallel alignments 
divided by the total resistance for the antiparallel alignment. 
As a function of voltage the variation of JMR is larger when 
one of the end metals is a paramagnet (cases c and d). For 
our plots we consider the ferromagnet FMi to be of iron, and 
FM2 is of cobalt. The parameters chosen are the same as in 
Fig. 2. 

modeli, we should mention that to see the spikes, the left 
or the right metal should be of a different ferromagnetic 
material. When all three metals are of the same ferro- 
magnetic materials, the steady state solution p"(t) is the 
same for parallel and antiparallel orientation of magne- 
tizations since the solution depends only on the ratio of 
the left and the right resistances. 

Comparing Figs. 4(a) - 4(b) with Figs. 4(c) - 4(d), 
we see that the variation of the JMR as a function of 
voltage is larger when one of the leads is non-magnetic. 
Following Eq. (|l|) and Eq. (||), which give the height 
of the plateau and the height of the first spike at zero 
temperature, we see that the height of the plateau for 
case (c) is lower than for case (a). On the other hand, 
the heights of the first spikes are same for both cases 
because the metals in the left and the center are same. 
Thus the variation of the JMR for case (c) is more than 
case (a). Cases (b) and (d) can be understood similarly. 

We plot the differential junction magnetoresistance as 
a function of the voltage in Fig. 5 for the same cases as 
in Fig. 4 and in the same parameter regime. Many of the 



4 



0.24 
0.20 
DC 0.28 

I 0-24 
0.00 
-0.40 



0.00 



' 1 ' 1 ' 1 




I: 





1 " 






I: 







n 



1 



CpV/e 



FIG. 5. Differential Junction magnetoresistance (DJMR) 
as a function of voltage for the same cases as in Fig. 4. We 
see spikes in the DJMR at those places where the current 
increases by a step. As in Fig. 4, the variation in the DJMR is 
larger as a function of voltage for cases (c) and (d) compared 
to cases (a) and (b). For case (c) the DJMR changes sign. 
Also compared to Fig. 4, the change in the DJMR is larger 
than the JMR. The metals are the same as used in Fig. 4 and 
the parameters are identical to Fig. 2. 



features of the DJMR are the same as for the JMR. The 
spikes are at the same places. The height of the spikes is 
also the same as in Eq. ( p^ . The DJMR at large voltages 
is obtained from the classical value with two resistors 
in series. Finally, the variation of the DJMR with the 
voltage is larger in cases 5(c) and 5(d) than in cases 5(a) 
and 5(b). There are, however, some differences. The 
DJMR can be negative, as shown in Fig. 5(c), while the 
JMR must be positive from Eq. (|l^) as long as i?^ < 
Rf^^ and Cn > Cl- Also, comparing Fig. (|) with Fig. 
(h), we see that the variation of DJMR as a function of 
voltage is larger than the variation of the JMR with the 
voltage. 



height of the spikes in general depends on the temper- 
ature, capacitances, and the resistances of the system 
whereas the height of the plateaus in the JMR vs. volt- 
age curves are given by the classical value obtained by 
adding two resistors in series. For asymmetric junctions 
where Cr > Cl , the spikes go up if the ratio of the resis- 
tances of the right junction for parallel and antiparallel 
alignment {R^/ R^) is greater than the same ratio of the 
resistances of the left junction (R^/R^^). Otherwise, the 
spikes go down. 

Many of the features of the DJMR are similar to the 
JMR: the spikes occur at the same places, the height of 
the first spike at zero temperature is the same, and the 
sign (up or down) of the spikes are the same. Also at high 
bias, the DJMR saturates to the classical value obtained 
by adding two resistors in series. However, there are few 
differences. The DJMR can be negative for some cases, 
while the JMR is positive as long as R^ < R^ and 
Cr > Cl- Finally, the variation of the DJMR is larger 
as a function of voltage than that of the JMR . 

In conclusion, double tunnel junction systems which 
have magnetic metals exhibit rich Coulomb blockade con- 
ductance vs. voltage curves. These curves provide a sig- 
nature of both the Coulomb charging effects and the spin 
polarization of the tunneling. By measuring the JMR 
or DJMR, one can extract information about both the 
capacitance charging energies and the spin polarization 
of the tunneling electrons, testing both theories of spin 
polarized tunneling and the Coulomb blockade. Further- 
more, there are cases where the JMR can be dramatically 
enhanced near one of the Coulomb blockade steps, mean- 
ing that there may be some applications of these effects. 
In any case, for small enough systems, both charging and 
spin-polarization effects will be important. 
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IV. CONCLUSION 



In this paper we have computed the junction magne- 
toresistance and the differential junction magnetoresis- 
tance for double tunnel junctions with magnetic metals in 
the Coulomb Blockade regime. As expected in a Coulomb 
Blockade problem, the steps in the current and the os- 
cillations in the total conductance are observed for both 
the parallel and antiparallel orientations of the magne- 
tizations, and the structure is reduced with increasing 
temperature. 

The JMR as a function of voltage shows spikes at the 
same places where the current increases by a step. The 
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